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ABSTRACT 


Aims. We aim to investigate the consequences of a fast massive stellar remnant — a black hole (BH) or a neutron star (NS) — encoun- 
tering a planetary system. 

Methods. We modelled a close encounter between the actual Solar System (SS) and a 2 Mo NS and a 10 Mo BH, using a few-body 
symplectic integrator. We used a range of impact parameters, orbital phases at the start of the simulation derived from the current SS 
orbital parameters, encounter velocities, and incidence angles relative to the plane of the SS. 

Results. We give the distribution of possible outcomes, such as when the SS remains bound, when it suffers a partial or complete 
disruption, and in which cases the intruder is able to capture one or more planets, yielding planetary systems around a BH or a NS. 
We also show examples of the long-term stability of the captured planetary systems. 
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1. Introduction 
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cg Observations have shown that the compact objects produced by 
— supernova (SN) explosions — neutron stars (NSs) or black holes 
(BHs) — can receive natal recoil impulses due to intrinsic asym- 
metries of their birth event. The best-fit observational velocity 
distribution of these kicks for NSs is a Maxwellian (Hansen & 
Phinney 1997; Jonker & Nelemans 2004) with a mean velocity 
pà about 300km s7! and a dispersion of about oy = 190km s™!. 
A BH could receive a kick of the same magnitude as a NS 
— (Jonker & Nelemans 2004; Repetto et al. 2017); although, more 
° recent investigations argue that lower velocities, for example 
= ~50kms~!, may be needed to explain the retention fraction of 
BHs in Galactic globular clusters (Peuten et al. 2016; Baumgardt 
& Sollima 2017; Pavlik et al. 2018), and that they are likely re- 
e -+ duced by the fallback of mass onto the BH remnant (e.g. Bel- 
. =, czynski et al. 2008; Fryer et al. 2012). 
In this work, we focus on the dynamical consequences of 
& a passage of a collapsed fast moving massive body close to 
© a planetary system using the current Solar System (SS) for re- 
alistic initial conditions. Such an event, which could more likely 
play out in the environment of a star cluster or an association, 
has hitherto not been fully examined. 

Published studies of the stability of planetary systems under- 
going stellar encounters have analysed isolated in-plane events 
involving low-mass single field stars or binaries (e.g. Li & 
Adams 2015; Malmberg et al. 2011) or have been sited in clus- 
ters with the impactor velocity taken from the virialised velocity 
dispersion using Monte Carlo or N-body procedures on fictitious 
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planetary systems (e.g. Cai et al. 2017; Wang et al. 2020). Our 
study differs in assuming (1) the real SS with all eight planets, 
(2) arbitrary off-plane interactions, (3) a massive intruder of high 
mass ratio, and (4) a broad range of impact parameters and en- 
counter velocities. The only study similar in spirit to ours, by 
Laughlin & Adams (2000), imposed an approximation of the SS 
and included high-velocity encounters with a low-mass binary 
system as the impactor. 


2. Simulation methods 


The set of nearly 2000 numerical simulations we performed in 
this study examined two encounter scenarios — BH and NS. We 
varied the impact parameter, by (defined in the plane of the SS), 
the encounter velocity, vy, the intruder mass, my, and incidence 
angle, 6; (relative to the invariable plane of the SS); the subscript 
‘T specifies the nature of the impactor: BH and NS. No relativis- 
tic effects were considered, so the designations BH and NS are 
formal and indicate a point mass. 

Our simulations were performed using REBOUND (Rein & 
Liu 2012), a symplectic N-body code, with a 15" order Gauss— 
Radau integrator IAS15 (Rein & Spiegel 2015), which was im- 
plemented for Python, and they were analysed using Python 
with NumPy (Harris et al. 2020) and Matplotlib (Hunter 2007). 
The system was initiated in isolation, ignoring the Galactic grav- 
itational potential and the very low stellar density in the Solar 
neighbourhood. The duration of the interaction is much lower 
than the Galactic orbital timescale of the whole system. The ini- 
tial conditions for specific encounter scenarios are described be- 
low in further detail. 


Article number, page 1 of 25 


A&A proofs: manuscript no. 40454-21 


2.1. Initial conditions 


We adopted the same range of impact parameters in both scenar- 
ios, 107! < by/pe < 10~° (the smallest, 0.3 au, is approximately 
the current orbit of Mercury). The SN kicks may not be the only 
reason for these encounters, they may also arise in a dynami- 
cal system (e.g. a star cluster or an association, as analysed by 
numerical models Spurzem et al. 2009; Hao et al. 2013; Zheng 
et al. 2015; Shara et al. 2016; Cai et al. 2017), so we chose ve- 
locities 10, 50, and 100kms7! in both scenarios. We also in- 
cluded 360 km s™! for the NS to account for possibly large ve- 
locity kicks or hyper-velocity stars. We also varied the impactor 
incidence angle, 6;, between 0° and 180° in 30° increments. In 
each setup, we launched the impactor from 0.5 pc towards the 
Sun and ended the simulation when the separation of these two 
bodies was again greater than 0.5 pc, since there was no effect at 
a larger distance 1 

The initial positions of the planets and the Sun were taken 
from NASA JPL HORIZONS system for Jan 1,2000 and then 
varied in each realisation by letting the SS evolve for a certain 
amount of time before the start of the simulation. We adopted 
nine pre-evolved states, evenly distributed from 0 to 165 years 
(approximately the orbital period of Neptune), producing an en- 
semble of possible encounters of the remnant with large planets. 
Such a coarse time step can initially create degenerate positions 
of the terrestrial planets but can be justified due to the following 
reasons: 


1. We let the massive body travel at different velocities from the 
same distance, so the positions of planets during the closest 
approach are different despite starting from the same pre- 
evolved state. 

2. We made several simulations without terrestrial planets. Us- 
ing deterministic symplectic methods allowed us to make 
one-to-one comparisons to the simulations with all planets 
and the same initial configuration. The results show that large 
gas planets are the dominant agents when it comes to the sta- 
bility of the SS or capture of planets. Low-mass terrestrial 
planets being retained, captured, or ejected is a secondary 
outcome that influences the orbital parameters of the final 
system should they be captured or retained. Hence, the ter- 
restrial planets behave essentially as test masses. We plot- 
ted one such comparison in Fig. 1 — the simulations are al- 
most identical except for the capture of Venus and Earth and 
a slightly different semi-major axis of Saturn. 

3. We also performed high-resolution calculations of the BH 
encounter, with a fixed incidence angle and a time step be- 
tween realisations of 0.5 year (i.e. 330 simulations in total, 
see Appendix A). 


2.2. Orbital parameters 


For analysing the simulation outputs, we distinguished planets 
that remained bound to the Sun from those captured by the in- 
truder. We counted all the orbiting bodies within a radius of 
200 au around the relevant central object that remained bound 
until the end of the simulation. The orbital parameters are de- 
fined as follows. The semi-major axis is given by 


-1 
2 v|? 
pa ed, (1) 
Irl GM 
' For two models, we also examined the longer-term stability of the 
captured systems, see Sect. 3. 
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where r and v are the radius and velocity vectors, respectively, 
in the co-moving frame of the central body (the Sun or the im- 
pactor), G is the Newtonian gravitational constant, and M is the 
mass of the central body. The orbital inclination is defined by 


i = arccos (n; / \hl) ; (2) 
where h = r x v in the co-moving frame, and the eccentricity 
vector is 

vxh r 


e= GM i ý 


where the orbital eccentricity is e = |e|. These values are plot- 
ted for each planet (labelled as Mercury ¥, Venus °, Earth @, 
Mars ©’, Jupiter 2, Saturn 5, Uranus 8, and Neptune ) that 
has been captured or retained in the top six rows of Figs. 2 & 3 
and in Figs. A.2, B.3-B.8, & C.3—-C.8. We emphasise that these 
plots show the distribution of orbital parameters, but not the ac- 
tual number of captured or retained planets, which is plotted, for 
example, in Figs. B.1, B.2, C.1, & C.2. 


3. Results 


Encounters with field stars have been simulated previously, but 
none specifically assumed that the intruder is a massive compact 
object. We therefore explored its interaction with a planetary sys- 
tem for arbitrary incidence angles which is likely in denser en- 
vironments, for example, star clusters (Varri et al. 2018) or the 
Galactic bulge, where thousands of NSs may be present (e.g. 
Fragione et al. 2018). The encounter velocities originating from 
SN kicks will be significantly higher than the virialised veloc- 
ities that have been used in published simulations (e.g. Wang 
et al. 2020). In fact, they are closer to what would be expected 
for a high-velocity star in the Solar neighbourhood (Laughlin & 
Adams 2000). An estimate of the gravitational capture distance, 
assuming that the intruder can be bound in the centre of mass 
frame, is given by Fgray1 = 2G (my + Mo) / v . This scales as 


Taav1 = 24x 1077 ( + all a 


-2 
-i , 4 
Mo a) ps 9 


where ‘I’ represents the impactor. For reference, the value for 
the lowest velocities in our simulations is Fgray,BH % 10-3 pe and 
Terav.NS © 2 X 10-4 pe. 

For bgus S Tgrav, We found alterations in both the incli- 
nation and eccentricity of virtually all retained planets regard- 
less of the incidence angle and encounter velocity (see Appen- 
dices B & C). The most efficient captures for both scenarios oc- 
curred in the range 60° <$ @gHNs S 120°, with up to eight planets 
for bgus S 1075 pe and the lowest impact velocity (see feap in 
Figs. 2 & 3). Eccentricities of the captured planetary orbits, ecap, 
were very high — around 0.8 for the BH and 0.65 for the NS. 

There was a difference between the prograde (Ogu,ns = 0°) 
and retrograde (@gHNns = 180°) encounters with the SS in both 
scenarios. A low-velocity NS showed no directional preferences 
for planetary capture, but it was most destructive for a retrograde 
fly-by with bys < 107% pc; in 60 % of the cases, the SS did not 
survive and in the remaining 40 % the Sun was left with only one 
terrestrial planet (see fret and fcap in Fig. 3). In contrast, a low- 
velocity BH was generally more disruptive and also more effi- 
cient in capturing planets for a prograde encounter (see fret and 
Seap in Fig. 2). 

The greatest difference between the NS and BH case was 
found in the change of the orbital inclinations of the retained 
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planets, iret, as the BH is more likely to completely flip the orbits 
of even massive planets. Based on the positions of the planets 
during the encounter, both intruders were able to capture planets 
on prograde, retrograde, or 90° orbits while approaching from 
OgH.Ns = 0° (see icap). It was, however, impossible for the NS to 
capture a retrograde planet arriving from Oys = 180°. 

Although the capture process is a many-body interaction, the 
final orbital angular momentum and semi-major axis of the cap- 
tured planet scales directly with the mass of the intruder. We also 
evaluated the three-body interaction of the Sun, Earth, and BH. 
These did not produce the outcomes we found in the complete 
simulations, as expected (see, e.g. Binney & Tremaine 2008; 
Heggie & Hut 2003). This three-body interaction is effectively 
a two-body encounter because of the large ratios of masses, and 
the ratio of the impact parameter and the semi-major axis of the 
Earth’s orbit. Consequently, the outcomes of our simulations re- 
quire the full interaction of the massive planets. 

At the smallest impact parameters and lowest velocities, the 
semi-major axes of the inner retained planets increased more 
than those of gas giants (see are in Figs. 2 & 3 and Figs. B.3 
& C.3). It should be noted, however, that the most frequent out- 
come was an ejection of the outer planets. Unlike previous stud- 
ies (Wang et al. 2020), we found no hot Jupiters among the re- 
tained systems even for the lowest impact parameters for any 
incident angle or velocity. We have also compared our NS in- 
plane encounter with the strong stellar fly-by scenario of Malm- 
berg et al. (2011). They used a schematic system containing only 
the giant planets, a 1.5 Mo intruder, 100 au for the impact param- 
eter, and 1 kms"! as the velocity at infinity. Our closest model is 
byns = 1074 pc and vys = 10km s™!. They found that x31 % of 
their systems ejected a planet immediately after the fly-by (listed 
in their Table 1) and we found a similar number, 25 %. In con- 
trast, they stated that ~28 % of their systems remained stable for 
up to 100 Myr, while we have no intact systems for any incidence 
angle. This difference is due to our higher impactor velocity and 
mass. Malmberg et al. (2011) do not provide the statistics on 
capture by the intruder, but we see 245 % cases where the NS 
has captured at least one planet, so the lack of stable systems is 
partly accounted for by this phenomenon rather than planetary 
ejection. 

We also performed two long-term (>1 Myr) integrations, one 
for each impactor, to examine the stability of the captured sys- 
tems. Both had captured Jupiter and three terrestrial planets. The 
NS was an orthogonal incidence, while the BH was from the 
high-resolution study described in Appendix A. In both cases, 
all planets were captured in highly eccentric orbits and remained 
bound to the impactor. The Jovian orbit was stationary while 
the terrestrial planets displayed the complex behaviour plotted in 
Fig. 4. We stress that these orbital histories are randomly chosen 
examples. Both cases showed long-term inclination oscillations 
between pro- and retrograde orbits, with the NS system being 
more regular and it also showed a larger excursion in eccentrici- 
ties than the BH case. In the former, Mercury nearly circularised. 
In contrast, the BH planetary system had a lower bound in the ec- 
centricity of all planets of ~0.5. In the NS case, the semi-major 
axes remained approximately constant for all planets while in the 
BH system, the terrestrial planets experienced abrupt jumps that 
never settled down. This may indicate the onset of longer-term 
chaotic behaviour. Further discussion is beyond the scope of this 
letter and a more extensive analysis will be a subject of a future 
study. 

Planetary capture by a NS is particularly interesting be- 
cause pulsars have been observed to posses planetary systems 
(even with multiple planets) — such as PSR 1257+12 (Wolszczan 


& Frail 1992) and PSR B1620—26 (Arzoumanian et al. 1996). 
While it is likely that these systems are relics (as, e.g. Wang et al. 
2006), our models provide an additional formation channel — the 
outcome of a close encounter of a NS with an ordinary planetary 
system. The mean rate for such encounters which could lead to 
a capture is about 3x10~7(n/pe~*) Myr™!, where n is the stellar 
density. For instance, this is 20.03 Myr! in star clusters which 
means acquiring a planet by capture could occur within the ac- 
tive lifetime of a pulsar ($108 years depending on the period). 
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Fig. 1. Sample realisation from the high-resolution study (see Appendix A) with vgy = 10kms~! and bgy = 107 pc; (a) with all planets and (b) 
without terrestrial planets. Both panels show the trajectories of the planets with coloured lines ( , Venus Ħ, Earth @, Mars ©', Jupiter 2, 
Saturn 5, Uranus 8, Neptune *¥), the Sun (yellow dots), and the BH (black dots). The left panels are plotted in the coordinates fixed to the initial 


position of the Sun (xo, yo, Zo), in which the impactor arrived from the right-hand side, and the right panels show a detailed view of the final system 
captured by the BH in its co-moving coordinates (gy, YBH, ZBH)- 
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Fig. 2. Outcomes of the BH scenario for vgy = 10kms™!. For each impact parameter (coded by colours and symbols) and incidence angle 
(columns), we plotted the final distribution of semi-major axes, eccentricities, and inclinations of planets that were either retained by the Sun or 
captured by the BH (top six rows, for each planet from Mercury ¥ to Neptune '?). The bottom two rows show the fraction of retained and captured 


planets. Each variation in the initial parameters was averaged over nine pre-evolved states of the SS, see Sect. 2. Detailed plots for all vgy and Og 
are in Figs. B.1-B.8. 
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Fig. 3. Outcomes of the NS scenario for vns = 10kms~!. For each impact parameter (coded by colours and symbols) and incidence angle 
(columns), we plotted the final distribution of semi-major axes, eccentricities, and inclinations of planets that were either retained by the Sun or 
captured by the NS (top six rows, for each planet from Mercury ¥ to Neptune '¥). The bottom two rows show the fraction of retained and captured 


planets. Each variation in the initial parameters was averaged over nine pre-evolved states of the SS, see Sect. 2. Detailed plots for all vys and 6ys 
are in Figs. C.1-C.8. 
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Fig. 4. Comparative long-term evolution of captured planetary systems. Left, NS: bys = 1076 pc, @6xs = 90°, vns = 50kms™!; Right, BH: 


bgy = 10° pc, 634 = —16°, vgy = 10 km s~}. The capture event is marked by a circle for each planet. The semi-major axis of Jupiter is not plotted 
in the top panels, but was 11.3 au in the NS system and 24.5 au in the BH system. See text for discussion. 
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Appendix A: High-resolution BH simulations 


Our study was initiated asking what would happen if a BH, ex- 
pected to form in the Betelgeuse SN, received a kick that would 
direct it towards the Solar System. Exploring this required spe- 
cific initial conditions linked to the progenitor and a high time 
resolution study. We found that the lower resolution generic 
cases that bounded these conditions were consistent. 

For what we hereafter call ‘the Betelgeuse case’, the remnant 
had mg = 10 Mo and was launched from the ecliptic latitude of 
the star, — 16.027° (J2000.0). We used this as the incidence angle 
in all of the high-resolution simulations. In the following discus- 
sion, we use the subscript « for ‘T to label all the parameters. We 
used encounter velocities 10kms=! < ve < 190kms~!. These 
were set up in the barycentric coordinates of the SS and, there- 
fore, are not the actual kick velocity of œ Ori which has to be 
corrected for the stellar radial velocity (Famaey et al. 2005), that 
is, the kick would be ~ 22 kms! higher. 

Our principal results are shown in Fig. A.2 and a sample re- 
alisation is shown in Fig. la. For by > 107° pc, the SS remained 
stable regardless of vg (see fret). The planetary orbits remained 
unchanged except for Mercury, whose inclination increased to 
5° (see dret, ret, Let). At these distances, the œ Ori remnant did 
not significantly alter the trajectory of the SS in the Galaxy — 
the maximum deviation of 0.8 km s7! (i.e. ~10 % of the peculiar 
velocity of the Sun) was obtained for ve = 10 kms! — with a de- 
creasing impact parameter, this effect increased (see Tab. A.1). 

For 10-2 > ba/pe 2 1073, the Sun still retained all planets if 
the remnant’s velocity was at least as great as the Earth’s orbital 
velocity. If it is very slow, vg = 10km s™!, it takes ~100 yr to 
traverse the SS, which is enough time to significantly alter the 
eccentricities of the two outermost planets. This either unbinds 
Neptune or results in its capture. 

For the highest velocity encounters, no planets were captured 
at any impact parameter above 107 pc. For all impact parame- 


Table A.1. Deviation (in absolute values) from the peculiar velocity of 
the Sun (Dehnen & Binney 1998), i.e. U = 10.00 kms"! (towards the 
Galactic centre), V ~ 5.25 km s™! (in the orbital motion of the Sun), and 
W = 7.17 kms™! (perpendicular to the Galactic plane). 


va /kms! > 190 100 50 29 10 
Ll ba/pe | kms"! 
107! AU 0.00 0.00 0.01 0.01 0.01 
AV 0.00 0.00 0.01 0.02 0.05 
AW 0.00 0.00 0.00 0.00 0.00 
10-2 AU 0.02 0.03 0.06 0.10 0.32 
AV 0.04 0.08 0.16 0.27 0.80 
AW 0.00 0.00 0.00 0.00 0.03 
1073 AU 0.15 0.29 061 1.08 4.77 
AV 0.42 0.81 1.61 2.66 5.78 
AW 0.00 0.00 0.02 0.10 2.20 
1074 AU 1.58 3.20 7.86 16.02 13.50 
AV 4.25 7.98 14.84 18.13 1.01 
AW 0.04 0.27 2.07 7.77 11.40 
1075 AU 19.23 47.47 63.19 40.76 13.00 
AV 40.43 57.68 17.15 4.04 4.28 
AW 3.83 21.87 46.66 34.55 11.90 
1076 AU 211.61 136.01 65.98 39.06 12.90 
AV 100.80 9.95 18.65 12.95 4.61 
AW 143.42 114.25 59.49 35.78 11.90 
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Fig. A.1. Frequency weighted capture cross section for the Betelgeuse 
scenario. The data points were calculated from our simulations. See text 
for details. 


ters by < 107° pc, the SS was considerably altered, especially 
regarding the number of planets retained and their orbital ec- 
centricities, as shown in Fig. A.2 (fret and eret). There were no 
discernible effects above impact parameters 107? pc, except for 
the 10 km s™! encounter velocity. 

The capture cross section resulting from the simulations in- 
cluding the entire system is Leap S 3x 107° pc”, which agrees 
with the simple estimate from Eq. (4). The number of captured 
planets depended sensitively on vg. At 10km s™!, all of the en- 
counters with by < 107° pe resulted in a capture of at least 
one planet in more than 15 % of the simulations. In contrast, 
at 50km s™!, only 1/330 simulations resulted in the capture of 
a single planet (Neptune, see Fig. A.2) with the maximum cross 
section Leap ¥ 3x1078 pc?. The number of planets captured was 
Np < 7, all of which had high orbital eccentricities and incli- 
nations. Regarding disruptions, at least 55 % of the simulations 
with by < 10-4 resulted in the ejection of at least one planet re- 
gardless of vg. We stress that these values are frequencies, not 
probabilities, because our models are not Monte Carlo simula- 
tions. As shown in Fig. A.1, we find a scaling relation for the 
a Ori capture cross section weighted by the frequency of models 
producing at least one capture of 
Scapa = 7X10 (va /kms“!) ” pe?. (A.1) 

As a final point regarding the Betelgeuse scenario, although 
a disruptive encounter of its remnant with the SS is highly un- 
likely, even a more probable distant encounter would not be par- 
ticularly pleasant. The models we present did not include any 
minor bodies (not even planetary moons). Their masses are neg- 
ligible relative to the planets or the Sun, so they have no effect 
on the dynamical stability of the SS, but it is certain that the 
close passage of a massive body will alter their positions and 
velocities. We performed a few simulations with Halley-type 
comets (a ~ 200 au), but even for an impact parameter as low 
as 107? pe, their orbits in the SS barycentric frame were only 
marginally altered. Nonetheless, a distant passage of a massive 
object, 0.1—1 pc, could strip away part of the Oort cloud and also 
produce an excess of comets with very long periods and high ec- 
centricities that could potentially cross Earth’s orbit, as found in 
previous studies of stellar fly-bys (Dybczynski 2006; Berski & 
Dybczynski 2016). 
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Fig. A.2. Outcomes of the Betelgeuse scenario. For each impact parameter (coded by colours and symbols) and encounter velocity (columns), we 
plotted the final distribution of semi-major axes, eccentricities, and inclinations of planets that were either retained by the Sun or captured by «œ Ori 


(top six rows, for each planet from Mercury ¥ to Neptune '?). The bottom two rows show the fraction of retained and captured planets. All data 
were averaged over 330 simulations. 
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Appendix B: Additional figures: BH simulations 


Fig. B.1. Capture of planets by a BH. Frequency of capture, fiap, of a given number of planets, Np, by the BH for the specified set of initial 
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conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, 0gp). 
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Fig. B.2. Planet retention after a BH encounter. Frequency of retention, fre, of a given number of planets, Np, by the Sun after the BH fly-by for 
the specified set of initial conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, Ogy). 
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Fig. B.3. Distribution of semi-major axes of planets captured by a BH. Each data point represents the mean value of the semi-major axis of 
a specific planet (Mercury ¥ to Neptune '}) if it was captured by a BH. The data were compiled from all pre-evolved states of the SS for a given set 
of initial conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, @gy) — i.e. the panels do not show how many planets 
were captured in total (see Fig. B.1 for reference). 
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Fig. B.4. Distribution of semi-major axes of planets retained by the Sun after a BH encounter. Each data point represents the mean value of the 
semi-major axis of a specific planet (Mercury ¥ to Neptune '}) if it was retained by the Sun. The data were compiled from all pre-evolved states 
of the SS for a given set of initial conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, 6g) — i.e. the panels do not 
show how many planets were retained in total (see Fig. B.2). The former SS is shown for reference. 
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Fig. B.5. Distribution of eccentricities of planets captured by a BH. Each data point represents the mean value of the eccentricity of a specific 
planet (Mercury ¥ to Neptune '}) if it was captured by a BH. The data were compiled from all pre-evolved states of the SS for a given set of initial 
conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, 6g) — i.e. the panels do not show how many planets were 
captured in total (see Fig. B.1 for reference). 
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Fig. B.6. Distribution of eccentricities of planets retained by the Sun after a BH encounter. Each data point represents the mean value of the 
eccentricity of a specific planet (Mercury ¥ to Neptune '}) if it was retained by the Sun. The data were compiled from all pre-evolved states of the 
SS for a given set of initial conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, Ogy) — i.e. the panels do not show 
how many planets were retained in total (see Fig. B.2). The former SS is shown for reference. 
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Fig. B.7. Distribution of inclinations of planets captured by a BH. Each data point represents the mean value of the inclination of a specific planet 
(Mercury ¥ to Neptune '}) if it was captured by a BH. The data were compiled from all pre-evolved states of the SS for a given set of initial 
conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, 6gy) — i.e. the panels do not show how many planets were 
captured in total (see Fig. B.1 for reference). 
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Fig. B.8. Distribution of inclinations of planets retained by the Sun after a BH encounter. Each data point represents the mean value of the 
inclination of a specific planet (Mercury ¥ to Neptune '}) if it was retained by the Sun. The data were compiled from all pre-evolved states of the 
SS for a given set of initial conditions (impact parameter, bgy, encounter velocity, vgy, and the incidence angle, Ogy) — i.e. the panels do not show 
how many planets were retained in total (see Fig. B.2). The former SS is shown for reference. 
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Appendix C: Additional figures: NS simulations 
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Fig. C.1. Capture of planets by a NS. Frequency of capture, fcap, of a given number of planets, N,, by the NS for the specified set of initial 
conditions (impact parameter, bys, encounter velocity, uns, and the incidence angle, yns). 
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Fig. C.2. Planet retention after a NS encounter. Frequency of retention, fre, of a given number of planets, Np, by the Sun after the NS fly-by for 
the specified set of initial conditions (impact parameter, bys, encounter velocity, vns, and the incidence angle, yns). 
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Fig. C.3. Distribution of semi-major axes of planets captured by a NS. Each data point represents the mean value of the semi-major axis of 
a specific planet (Mercury ¥ to Neptune '¥) if it was captured by a NS. The data were compiled from all pre-evolved states of the SS for a given set 
of initial conditions (impact parameter, bys, encounter velocity, vns, and the incidence angle, ns) — i.e. the panels do not show how many planets 
were captured in total (see Fig. C.1 for reference). No planets were captured in the fastest scenario. 
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Fig. C.4. Distribution of semi-major axes of planets retained by the Sun after a NS encounter. Each data point represents the mean value of the 
semi-major axis of a specific planet (Mercury ¥ to Neptune '}) if it was retained by the Sun. The data were compiled from all pre-evolved states 


of the SS for a given set of initial conditions (impact parameter, bys, encounter velocity, vns, and the incidence angle, yns) — i.e. the panels do not 
show how many planets were retained in total (see Fig. C.2). The former SS is shown for reference. 


Article number, page 21 of 25 


A&A proofs: manuscript no. 40454-21 


bys/pe: © 1078 A 10-4 v 10 o 10-6 


Uys = 360kms7! vys = 100 km s7} vys = 50km s7! vys = 10km s7} 
1.00 y z z= 
v 
= o a. B ï g aH 
5 Oo a = Oo foo 
0.754 4 4 a [m] v 
o 2 
2 0.504 J l 
vo n 
2 
D 
o 
S 
n a) 
5 ll 
v n 
2 
D 
2 
7 5 
s | 
w n 
zZ 
D 
o 
= 
a Ss 
3 l 
vo n 
Z 
D 
aS 
N 
a kan 
< ll 
Zz 
D 
aS 
Yo) 
2 Smal 
E I 
A 
D 
o NOgy a v 
0.754 4 J ʻo Ea vY oyða 
v A |© 
a re 
E] 
= 0.504 : | 1 
n 
A |e 
0.25 4 4 4 
0.00 


>PatTr*4H EYP FRAGTAHEYP FRAGTAHEYP FRAGTASKEP 

Fig. C.5. Distribution of eccentricities of planets captured by a NS. Each data point represents the mean value of the eccentricity of a specific 
planet (Mercury ¥ to Neptune '?) if it was captured by a NS. The data were compiled from all pre-evolved states of the SS for a given set of initial 
conditions (impact parameter, bys, encounter velocity, vys, and the incidence angle, ns) — i.e. the panels do not show how many planets were 
captured in total (see Fig. C.1 for reference). No planets were captured in the fastest scenario. 
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Fig. C.6. Distribution of eccentricities of planets retained by the Sun after a NS encounter. Each data point represents the mean value of the 
eccentricity of a specific planet (Mercury ¥ to Neptune 'P) if it was retained by the Sun. The data were compiled from all pre-evolved states of the 
SS for a given set of initial conditions (impact parameter, bys, encounter velocity, vys, and the incidence angle, Ons) — i.e. the panels do not show 
how many planets were retained in total (see Fig. C.2). The former SS is shown for reference. 
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Fig. C.7. Distribution of inclinations of planets captured by a NS. Each data point represents the mean value of the inclination of a specific planet 
(Mercury ¥ to Neptune '}) if it was captured by a NS. The data were compiled from all pre-evolved states of the SS for a given set of initial 
conditions (impact parameter, bys, encounter velocity, vys, and the incidence angle, Oys) — i.e. the panels do not show how many planets were 
captured in total (see Fig. C.1 for reference). No planets were captured in the fastest scenario. 
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Fig. C.8. Distribution of inclinations of planets retained by the Sun after a NS encounter. Each data point represents the mean value of the 
inclination of a specific planet (Mercury ¥ to Neptune '}) if it was retained by the Sun. The data were compiled from all pre-evolved states of the 
SS for a given set of initial conditions (impact parameter, bys, encounter velocity, vys, and the incidence angle, Ons) — i.e. the panels do not show 
how many planets were retained in total (see Fig. C.2). The former SS is shown for reference. 
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